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ABSTRACT 



Aims. The aim of this work is to characterize the average Fe K emission properties of active galactic nuclei (AGNs) in the source 
rest-frame. 

Methods. We selected a sample of 248 AGNs from the 2XMM catalog, covering a wide redshift range < z < 5 and with the EPIC- 
PN net 2-10 keV rest-frame counts >200 and power law photon indices in the range 1.5-2.2. We employed two fully independent 
rest-frame stacking procedures to compute the mean Fe K profile. The counting statistics for the integrated spectrum is comparable to 
the one available for the best studied local Seyferts. To identify the artifacts possibly introduced by the stacking procedure, we have 
carried out simulations. 

Results. We report that the average Fe K line profile in our sample is best represented by a combination of a narrow and a broad line. 
The equivalent widths of the narrow and broad (parametrized with a diskline) components are ~30 eV and ~100 eV, respectively. 
We also discuss the results of more complex fits and the implications of the adopted continuum modeling on the broad line parameters 
and its detection significance. 

Key words, galaxies: active - X-rays: galaxies - quasars: emission lines 



1. Introduction 

Fluorescent Fe Kff lines in active galactic nuclei (AGNs) are a 
potentially unique tool to probe the innermost regions around a 
black hole. The line properties, such as the peak energy, intensity 
and profile carry important diagnostic information about the dy- 
namics and physics of the region where the emission originates. 
The measured line energy can be used to infer the ionization state 
of the line emitting matter, whereas the line equi valent width 
(EW) indicates the amount of fluorescing material (iFabian et al.l 
12000) 



The last decade has witnessed a significant improvement in 
our knowledge of the narrow Fe Ka line properties in AGNs. 
Based on data collected with high-sensitivity X-ray satellites 
XMM-Newton, Chandra and Suzaku the ubiquitous presence 
of a neutr al, narrow Fe Ko- line in AGNs spectra has been es- 
tablished ( Yaqoob & Padmanabhan 2004;' Guainazzi et al.ll2006l: 
Reeves et al.. ,2006; .Nandra et al.. .2007.) . The behavior of the 
Fe Ka line as a function of the X-luminosity and redshift has 
also been examined. Compelling evidence of the inverse cor- 
relation between the neutral, narrow Fe Ka line EW and X- 
ray luminosity, known as the X-ray Bal dwin effect (Iwasawa- 
Taniguci effect), has been reported (|I\yasawa & Taniguchil 



The line profile is determined by several physical parame- 
ters. The Fe K line in AGNs is assumed to be produced through 
X-ray irradiation of optically thick matter, s uch as the molec 
ular torus and/or t he ac c retion disk (,Guilbert & Rees 19881 
Lightma n & White! Il988t iGeorge & Fabiar3 Il99ll: iMattetal 



19931; iNandra et al.lll997l: iPage et al.ll2q04l: iBianch i et al.ll2007 
'Chaudha rvetal.ll2010t IShu et al.l 1201 oT Ho wever, the EW is 



1991.) . If the Fe Ka line arises in distant material like the molecu- 



lar torus e nvisaged in orien tation-dependent unification scheme 
for AGNs (I Antonuccil 1 9931) . the resulting emission line profile is 
narrow. In contrast, the Fe Ka line originating in the inner accre- 
tion disk is distorte d by Doppler and gravitational e ffects and be- 
comes asymmetric (IFabian et al.lll989HLaol^ll99Tl) . Hence, rela- 
tivistic lines offer a robust way to measure accretion disk prop- 
erties, such as the radial extent, the emissivity profile, the incli- 
nation angle of the disk to the observer's line of sight and the 
black hole spin (for the reviews see [Reynolds & Nowak 20031 
IFabian & Miniuttill2005t lMilleill2007l:lTurner & Milleni2009i) . 



observed to be indepe ndent of the redshift (iBrusa et al ] l200l 
IChaudharv e"tai]|2010t) . 

Since the first unambiguous detection of relativistically 
broadened Fe Ka emission in the X-ray s pectrum of MCG-6 - 
30-15 observed with the ASCA satellite ffanaka et anil995l) . 
sincere efforts have been made to test for the presence of broad 
Fe Ka emission lines in AGNs. Recent X-ray observatories like 
Chandra, XMM-Newton and Suzaku have revealed that statis- 
tically significant broad Fe lines can be detected in relatively 
deep observations. Using a sample of 102 AGNs o bserved with 
XMM-Newton at z < 0.5, iGuainazzi et al.l (l2006l) detected rel- 
ativistic lines in about 25% of the sample objects (^50% when 
the "well-exposed" spectra with >10 000 net counts in the 2-10 
keV band are considered). They also stacked residuals in four 
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equally populated luminosity classes and find that the average 
EWs are <150 eV for all luminosity classes. iNandra et aP (l2007h 
performed a spectral analysis of a sample of 26 nearby Seyfert 
galaxies (z < 0.05) observed by XMM-Newton with a minimum 
of 30 000 net counts in the EPIC-PN spectrum in the 2-10 keV 
band. They found evidence of broad line emission in around 65% 
of the sample with a typical EW of ~77 + 16 eV, w hen fitted with 
a broad Gaussian. Ide La Calle Perez et alj ( 1201 Ol) carried out a 
systematic and uniform analysis of 149 local radio quiet Type 
1 AGNs and found strong evidence of a relativistic Fe Ko' line 
in 36% of the sour ces, with an average EW of 100 eV (see also 
IPatrick etalj|2010l) . 

The X-ray spectra of the distant AGNs are often photon 
starved. The low photon counting statistics in such AGNs spec- 
tra prevents an appropriate modeling of the underlying contin- 
uum and therefore the significance of the broad line component 
and its measured parameters suffer from considerable uncertain- 
ties. The stacking analysis, a powerful technique to improve the 
counting statistics, can be applied to gain an insight into the av- 
erage properties of faint AGNs. 

Using a large sample of 3 52 AGNs de t ected in the Chandra 
Deep Field North and South, iBrusa et aP (l2005h computed the 
average spectra by stacking X-ray counts in the observed frame 
in seven red shift bins over the rang e 0.5 < z < 4. In a pre- 
vious paper dChaudharv et al.|[2010l hereafter CHIO), we as- 
sembled a sample of 507 sources from the 2XMM catalog at 
high galactic latitude {\BII\ > 25 degrees), with the sum of the 
EPIC-PN and EPIC-MOS 0.2-12 keV counts greater than 1000 
and covering the redshift range < z < 5. We stacked X-ray 
spectra in the observed frame in narrow redshift bins and inves- 
tigated the fundamental properties of the integrated spectrum, 
such as the ubiquity of the Fe Ko- emission in AGNs and the 
dependence of the spectral parameters on the X-ray luminosity 
and redshift. However, the analyses of iBrusa et al.l (l2005h and 
CHIO were limited to the narrow line properties as the broad 
line parametrization was hampered by redshift smearing in the 
ob served frame stacked spect ra. 

IStreblvanska et al] (l2005h constructed an average rest-frame 
spectrum of 104 X-ray sources (53 Type 1 and 41 Type 2 AGNs 
and 10 galaxies) covering < z < 4.5 in redshift space and 
reported a clear relativistic line in the average rest-frame spec- 
tra with an EW of ~560 e V and ~460 eV for the Type 1 and 
Type 2 AGNs, respectively. iLonginotti et al.l (l2008l) stacked X- 
ray spectra from a local sample of 157 AGNs and found that the 
EW of the broad relativistic line is never higher than 100 eV, 
either when stackin g the whole sample or different sub-samples. 
ICorral et"aD (l2008h computed the mean Fe emission from a large 
sample of more than 600 Type 1 AGNs spanning a redshift range 
up to ~3.5. They detected significant narrow Fe Ka emission line 
around 6.4 keV with an EW~90 eV and found no compelling 
evidence of any significant broad relativistic emission line in the 
average spectrum (see also lMao et al 

While the presence of a narrow line component is confirmed 
by all independent analyses, the evidence of the presence of a 
broad line is more debated. The complex nature of AGNs spec- 
tra, and, in particular, the degeneracy between the spectral pa- 
rameters and the dependence of the line EW and shape over the 
fitted (absorbed) continuum, hampered so far a comprehensive 
interpretation of the physics behind the observed features. As an 
example, it has also been proposed that the observed red-wing 
can also be ascribed to complex (clumpy) absorption at least for 
MCG-6-30-15 (Miller et al 2009). It is therefore critical to have 
multiple independent procedures to determine the Fe line param- 
eters. 



In this paper, we study a small but well-defined sample of 
AGNs taken from CHIO, and concentrate on the mean proper- 
ties of the Fe K line in the rest-frame using two fully indepen- 
dent procedures. The paper is organized as follows: Section|2]de- 
scribes the selection criteria and properties of the sample. Sect. 
[3] outlines the stacking procedures. The results of the spectral 
analysis are presented in Sect. |4] and are discussed in Sect. |5] 
Throughout this work a cosmology with Q„, - 0.27, Q.^ = 0.73 
and Hq - 70 km s"' Mpc"' is used. Spectral analysis was car- 
ried out with XSPEC (version 12.6.0); errors are reported at the 
68% confidence level for one interesting parameter (Ax~ = 1.0). 
All energies refer to the rest-frame unless otherwise specified. 



2. The sample 

We start with the CHIO sample of 507 AGNs at high galactic 
latitude (\BII\ > 25 degrees), with the sum of the EPIC-PN and 
EPIC-MOS 0.2-12 keV counts greater than 1000 and covering 
the redshift range < z < 5. This sample was selected from the 
2XMM catalog, the second comprehensive catalog of serendip- 
itous X-ray sources from the European Space Agency's (ESA) 
XMM-Newton observatory. The catalog contains ~250 000 X- 
ray sour ce detections which relate to ~200 000 unique X-ray 
sources (IWatson et al.ll2009h . Redshifts were obtained from the 
Nasa's Extragalactic Database. For each object, source specific 
products such as the source and background spectra, images in 
the different energy bands and response files were retrieved from 
the XMM-Newton Science Archive. We refer the reader to CHIO 
for all the details on the selection of the sample and safety checks 
performed to assess the quality of the archival products. 

An essential requirement of the conventional rest-frame 
stacking analysis is the spectrum unfolding by appropriately 
modeling the underlying continuum. Therefore, to carry out a 
reliable spectral analysis of each spectrum, we refined the sam- 
ple of CHIO by considering only the AGNs with the EPIC-PN 
net 2-10 keV rest-frame counts >200 (-73% of the parent CHIO 
sample). As our main interest is in the Fe K region of the spec- 
trum, we performed a detailed spectral analysis of the individual 
spectrum in the rest-frame 2-10 keV energy range, excluding the 
5.5-7.0 keV Fe bancfl using XSPEC (ver. 12.6.0; Arnaud 199^. 
We fitted the data with an absorbed power law0. Only the EPIC- 
PN data has been used in this work. 

From the original CHIO sample of 507 sources, we further 
selected a sample of 248 AGNs having their power law photon 
indices as derived from the above mentioned fits in the range 1 .5 
< r < 2.2. The imposed threshold on F ensured that our results 
are stable against the continuum normalization and the contri- 
bution of the highly obscured AGNs with prominent Fe lines is 
minimum. The top row in Fig. [1] shows the distributions of the 
EPIC-PN net counts (left) and power law photon indices in the 
rest-frame 2-10 keV band (right) of the 248 sources used in this 
work (green cross-hatched) compared to the CHIO parent sam- 
ple (red solid line). The bottom row shows overlaid histograms 
of redshifts (left) and 2-10 keV luminosities, not corrected for 
absorption (right) for these two samples. We note that the 248 



' Spectral fits excluding 5-7 keV and 6-7 keV have also been per- 
formed. This did not change the results. 

^ We also tested a model with a ne utral reflection component 
(pexrav in XSPEC, see iMagdziarz & Zdzi arski 199^). However, our 
choice of the spectral energy range (2-10 keV) and the limited statistics 
in the majority of spectra (see Fig.[T]for sample properties) restrained a 
full parametrization of the reflected continuum. We therefore adopted a 
simple absorbed power law continuum for all sources. 
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Fig. 1. Top: distributions of the EPIC-PN net counts (left) and 2-10 keV power law photon indices (right) of the 248 sources selected 
for our study (green cross-hatched) compared with the CHIO parent sample (red solid line). Bottom: redshift (left) and rest-frame 
2-10 keV luminosity, no absorption correction (right) distributions of these two samples. 



sources subsample can be considered a fair representative of the 
CHIO sample in terms of the X-ray luminosity and redshift. 

3. Rest-frame stacking procedures 

3.1. Stacked ratio 

A commonly employed method to study the average Fe line 
properties is to derive the average of the ratio between the data 
and a simple continuum model for each s ource under study (e.g. 
iNandra et al.lll997l:lGuainazzi et al.ll2006h . However, several key 
factors should be considered in the stacking analysis when the 
sample sources have a wide redshift distribution. These include 
different redshift of each source and thus a different observed 
line energy, the energy dependent detector response is also dif- 
ferent for each source. 

Taking into account all these factors we adopted the follow- 
ing procedure for constructing an average ratio profile for our 
AGNs sample. 

• For each source, we determined the observed frame energy 
range corresponding to the rest-frame 2-10 keV band and 
the associated ungrouped channel information from the re- 
spective response files. Each spectrum, in the rest-frame 2- 
10 keV band, was then grouped in predefined bins of widths 
equal to 0.25/(1 + z) up to 8 keV and 1/(1 H- z) in the 8-10 
keV using the FTOOLS routine grppha, which corresponds to 



a rest-frame samphng of 0.25 keV and 1 keV in the 2-8 keV 
and 8-10 keV, respectively. Thus, we had 26 energy bins in 
each spectrum. 

• Given the limiting counting statistics and the choice of pre- 
defined binning we could not assure a minimum of 20 counts 
per bin. We therefore analyzed the binned, background- 
subtracted, 2-10 keV rest-frame spectra with an absorbed 
power law using the Cash statistic (Cash 1973) implemented 
in XSPEC. The best fit parameters (Nh and F etc. with errors) 
and ratio with respect to the best fit continuum model were 
saved. 

• These ratios were then summed and averaged for all the 
sources in each rest-frame energy bin. First, we derived the 
mean and standard deviation of the ratios in each energy bin. 
We then removed in each energy bin all the ratios deviating 
more than 3 times the standard deviation from the average 
value (3-sigma clipping). Figure |2] shows the resulting aver- 
aged ratio (red) comprising >95% of ratios in each energy 
bin compared with the mean ratio for the whole sample of 
the 248 sources (blue). We find that the two ratio profiles are 
consistent within the errors as expected. We note that the 3- 
sigma clipped averaged ratio (red) is comparatively lower at 
high energies (> 7 keV) as the contribution of the ratios de- 
viating significantly from the average value in these bins is 
removed. 
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Fig. 2. Comparison of the mean ratio profiles: the whole sam- 
ple comprising 248 sources (blue) and after applying 3-sigma 
clipping on the ratios in each energy bin (red). 
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Fig. 3. Ratio of the "averaged X-ray spectrum" with respect to a 
power law (green) overlaid with the mean ratio profile of the 248 
sources created using 3-sigma clipping (red). 



• The averaged ratio cannot be fitted directly in XSPEC and 
it must be converted to a flux spectrum taking into account 
the average underlying continuum, we therefore converted 
the averaged ratio created by applying 3-sigma clipping (red 
points in Fig. |2]i to a flux spectrum by multiplying the av- 
eraged ratio in each energy bin by the E factor, where E 
is the central rest-frame energy of the bin and r(=1.8) is the 
mean slope of our sample. We refer to this flux spectrum as 
the "averaged ratio flux spectrum". 

3.2. Stacked spectrum 

We followed the averaging procedure of Iwasawa et al. (2011) 
for constructing an "averaged X-ray spectrum". The advantage 
of this procedure is that it is independent of the continuum mod- 
eling and statistics used in the spectral analysis. In this procedure 
original data, in units of counts s"' keV ' , of each spectrum were 
saved after loading the spectrum grouped as described above. We 
also created the ascii tables of the ancillary response matrices 
(.arf) with the same binning as of the spectrum. The data were 
then divided by the average effective area in each energy bin, and 
thereafter normalized to the 3-5 keV continuum to ensure that 
the average spectrum is not dominated by the brightest objects. 
The final "averaged X-ray spectrum" is constructed by summing 
the normalized counts of all the different sources in each rest- 
frame energy bin and averaging for the total number of sources. 
A special effort is made to evaluate calibration issues and uncer- 
tainties (see Iwasawa et al. (201 1) for further details). The total 
number of the 2-10 keV counts for the 248 sources is ~ 198000. 

The ratio of the "averaged X-ray spectrum" with respect to a 
power law (green) over-plotted with the mean ratio profile cre- 
ated using 3-sigma clipping (red) is shown in Fig. [3] We notice 
that the two ratio profiles are fully consistent within the statis- 
tical errors. Spectral fitting of the "averaged X-ray spectrum" is 
addressed in more detail in Sect. 14.21 

3.3. Simulations 

To identify the artifacts possibly introduced by the stacking 
method, we carried out simulations using XSPEC. First, we cre- 
ated 100 simulated "ungrouped" line less spectra for each of the 
248 sources using their best spectral fit parameters (e.g, Nh, F 
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Fig. 4. Averaged ratio of the 248 spectra after applying 3-sigma 
clipping (red), mean simulated continuum (black) along with its 
Icr (green dotted line) and 3cr (blue solid line) confidence limits. 
The solid vertical line is drawn at 6.4 keV. 



and power law normalization etc.). Identical response and effec- 
tive area files were used for the real and simulated spectra. 

We then applied the same procedure detailed in Sect. l3.1l to 
these simulated spectra and computed the mean ratio profile of 
the 248 simulated spectra for each realization. Thus, we have 
100 realizations of the underlying continuum represented by the 
mean ratio of each realization. We then computed the distribu- 
tion of these simulated continua. The averaged ratio profile of 
the 248 spectra (red), mean simulated continuum (black) and its 
Icr (green dotted line) and 3cr (blue solid line) deviations in each 
energy bin are shown in Fig.|4] 

We note that some apparent deviations (e.g. upward turn at 
energies below 2.5 keV) present in the averaged ratio profile is 
also evident in the averaged simulated continuum. This in turn, 
confirms that simulations can identify very small systematic ef- 
fects of the order of 5% and the significance of any spectral fea- 
ture can be deduced with a reliable accuracy. We have investi- 
gated the possible reasons for this curvature at energies below 
2.5 keV. This effect can be attributed to a combination of the 
boundary effect (due to our spectral energy range selection) and 
counting statistics. Indeed, when the averaged ratio profile and 
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the averaged simulated continuum are computed for the sources 
with net counts >500, we do not find any curvature. 

4. Results 

4.1. "Averaged ratio flux spectrum" fitting 

We performed a detailed spectral analysis of the "averaged ratio 
flux spectrum" constructed from the averaged ratio profile as de- 
scribed in Sect. 13. II All the fit results are shown in the Appendix. 
We used the FTOOLS routine f lx2xsp to convert the ascii file of 
the "averaged ratio flux spectrum" to the fits format. To assess 
the improvement in the spectral fitting on including an additional 
spectral feature in the fit, we adopt the Ax^ criterion (e.g. Ax^ 
= 2.7 corresponds to 90% significance for the addition of one 
interesting parameter and Ax^ = 4.6 corresponds to 90% signifi- 
cance for two additional interesting parameters). Throughout the 
spectral analysis, errors and upper limits are quoted at the 68% 
confidence level for one parameter of interest. The errors in the 
Fe line EW are calculated using XSPEC eqw with err option. 

We first carried out the spectral analysis in the 2-10 keV en- 
ergy range using various models consisting of simple/complex 
continuum and one/two Fe Ka lines. However, in all the at- 
tempted fits, we noticed the presence of significant residuals of 
the order of 5% in the 2-3 keV band that can be ascribed to the 
systematic eff'ect identified by our simulations discussed in sec- 
tion l3.3l Moreover, the 2-3 keV excess resulted in high;^'^ values 
(e.g. x^ > 2) in most of the models explored. Therefore, to ob- 
tain more reliable spectral fits, we restricted the spectral analysis 
of the averaged ratio flux spectrum in the 3-10 keV band. 

We then concentrated on the characterization of the narrow 
Fe K line. To parametrize the narrow core we used a power law 
plus a narrow Gaussian line with the energy and width fixed 
at 6.4 keV and 0.01 keV, respectively. The line normalization 
was allowed to vary. This line model yielded an acceptable fit 
(x^/dof = 31.3/19) and the narrow line is detected at >99.995% 
significance level (Ax^ ~ 30 for one additional parameter), with 
an EW of 45 + 13 eV. We note that the residuals show excess 
emission of the order of a few per cent in the 6-7 keV energy 
range. 

To search for the presence of a broad component we replaced 
the narrow feature with a broad Gaussian by leaving the line 
width free to vary. The fit was considerably better than a sin- 
gle narrow Gaussian emission line model (x^/dof = 25.6/18), 
with resulting values of the broad line EW of 75 + 40 eV and a 
measured width of o" = 0.15;^q q4 keV. However, we see excess 
emission around 6.7 keV possibly indicating that the model with 
a single line component does not fully describe the data. This hy- 
pothesis is further supported by the observed improvement in the 
goodness of fit as a single Fe K line with free width is partially 
modeling the broad feature. 

We next fitted the data with a simple power law and a phys- 
ically motivated broad line with an emission profile from an 
accretion disk aro und a Schwarzschild black hole (diskline, 
iFabian et al.lll989h . We fixed several diskline parameters; the 
peak energy, Eoisk, at 6.4 keV, the inner disk radius, 7?,,,, outer 
disk radius, Rom, at 6 and 100 Rg, respectively, where Rg is the 
gravitational radius and the emissivity index, yS, at -2.5. The free 
parameters were the disk inclination, /, and normalization. The 
fit was statistically worse (x^/dof - 41.6/18) as compared to 
the previous ones, which can be attributed to the presence of 
significant residuals peaking at 6.4 keV. The errors in the fit pa- 
rameters except the diskline EW (128 + 41 eV) could not be 
inferred due to the poor fit. Thus, we concluded that a line model 





CO 




O 








X 




00 








en 

1 


cm- 


o 

X 




CD 


C/3 






CO 


x: 


Q. 


O 








X 


a? 




_c 




15 


CO 




o 


b 


X 







_l I I L 




J I I l_ 



_l I l_ 



10-3 2x10-3 
I (Narrow) (ph s-^ cm-^) 



3x10-3 



Fig. 5. Contour plots (68%, 90% and 99%) of the diskline ver- 
sus narrow component intensities for the "averaged ratio flux 
spectrum" . A simple power law is used to model the underlying 
continuum. The cross indicates the best fit values. 



consisting of two components should be used to model the Fe K 
band in our data. 

We therefore tested a two component line model consist- 
ing of two Gaussians to account for the narrow core and the 
broad component. The fixed parameters for the narrow com- 
ponent were the same as those in the single narrow Gaussian 
line fit mentioned above. The line energy of the broad feature 
was fixed at 6.4 keV, while the broad line width and normaliza- 
tions of both the lines were free parameters. With an improve- 
ment in goodness of fit (x^/dof = 22.8/17), the fit with two line 
components gives narrow and broad line EWs of 28 + 20 eV 
and 73 ± 43 eV, respectively. Using a simple power law contin- 
uum parametrization, a broad Gaussian Fe Ko' line is detected at 
~95% significance {Ax~ ~ 9 for two additional parameters), with 
the width of cr = 0.40;^^ keV. On modeling the broad feature 
with a diskline with fiit parameters of a single diskline fit, 
we achieved a similar quality fit (x^/dof = 22.7/17). The mea- 
sured EWs of the narrow and diskline components are 39+14 eV 
and 102^^j eV, respectively, while the disk inclination is 44. 6;^^ " 
degrees. Figure |5] plots contours (at 68%, 90% and 99% confi- 
dence interval) of the diskline versus narrow component in- 
tensities together with the best fit values derived for this fit (mo 
pow+gaussian+diskline) to the "averaged ratio flux spec- 
trum". We also checked the variation in the diskline EW by 
fixing the outer disk radius, to 1000 Rg as well to the maxi- 
mum value of 10^ Rg allowed by the diskline model. The mea- 
sured diskline EW is found to be within the 2 eV difference 
of the previous value. 

Given the robust detection of the narrow Fe Ka line in the 
"averaged ratio flux spectrum", a signature of reprocessing in 
cold, distant material, we a dded a neutral reflection comp onent 
using the pexrav in XSPEC dMagdziarz &Zdziarskil l 19951) . The 
pexrav model computes the spectrum resulting from an expo- 
nentially cut-off power law continuum incident on a slab of op- 
tically thick, neutral material, accounting for the Compton re- 
flection and bound-free absorption. However, the pexrav does 
not include any fluorescence emission lines. We use the pexrav 
to parametrize the reflected continuum only, as the illuminat- 
ing primary continuum is separately modeled. In this model, the 
pexrav continuum shape was tied to the mean photon index of 
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Fig. 6. Contour plots (68%, 90% and 99%) of the diskline ver- 
sus narrow component intensities for the "averaged X-ray spec- 
trum". A simple power law is used to model the underlying con- 
tinuum. The cross indicates the best fit values. 



the illuminating continuum, F (fixed at 1.8), abundances were 
assumed to be Solar, the inclination of the reflector to the ob- 
server's line of sight, cosi, and the cut-ofF energy, foldE, were 
fixed at 0.90 and 200 keV, respectively. The solid angle of the re- 
flector Q/27r was fixed to 1. The reflection fraction of the distant 
neutral reflector, Roist, is derived from the ratio of the reflection 
component normaUzation to the power law normalization. The 
errors in the Rms t are propagated using the standard error prop- 
agation formula ( Bevingtonll 19691) . To account for both the line 
features we first used two Gaussian lines with parameters set as 
in the power law continuum modeling. 

Overall, the fit statistic is good (x^/dof = 21.6/17) with/?/),,, 
= 0.27 ±0.17 and the naiTowline EW of 28+20 eV. With this con- 
tinuum parametrization, the detection significance of the broad 
Gaussian line with cr - 0.37;^g 2Q is marginal {Ax^ ~ 4 for two 
additional parameters). We therefore inferred an upper limit of 
105 eV for the broad line EW. When the broad component is 
modeled by a diskline profile, the fit statistic is similar to the 
previous fit, with (x^/dof = 21.5/17) and the EWs of the nar- 
row and diskline are 37 + 14 eV and <137 eV, respectively. The 
fit results are summarized in Table [1] The observed variation in 
the broad line detection significance as a function of the adopted 
continuum parametrization is noticeable. 

4.2. "Averaged X-ray spectrum" fitting 

We have also carried out a detailed spectral analysis of the "av- 
eraged X-ray spectrum" in the 3-10 keV energy range using 
the four models discussed above with 2 line components. We 
mention that in the "averaged X-ray spectrum" fitting, the free 
and fixed parameters of the two lines (parametrized by either 
Gaussians or a narrow and a diskline) are set in accordance 
with those in the "averaged ratio flux spectrum" fitting. First, 
the underlying continuum was modeled with a power law and 
two Gaussian lines representing the narrow and broad features 
were included. The fit statistic was reasonably good (x^/dof - 
20.6/17) with the power law photon index, F, of 1.64 + 0.03 
and narrow line EW of 28 + 10 eV, while the broad line EW 
is 92 + 35 eV. The broad Gaussian line is detected at the ~99% 
confidence level {Ax^ -II for two additional parameters), being 
in an excellent agreement to the corresponding fit of the "aver- 



aged ratio flux spectrum". When the broad feature is modeled 
by a diskline, the fit statistic is exceflent (x^/dof = 17.4/17) 
and we recovered the narrow line EW = 37 + 8 eV, diskline EW 
= 128 + 36 eV and the disk inclination ; = 44.4 + 2.4 degrees. 
The diskline is detected at a robust significance level [Ax^ ~ 
14 for two additional parameters). Figure |6] shows contours (at 
68%, 90% and 99% confidence interval) of the diskline ver- 
sus narrow component intensities together with the best fit val- 
ues derived for this fit (mo pow+gaussian+diskline) to the 
"averaged X-ray spectrum". 

The additional spectral fits with a simple power law and the 
reflection component (modeled with pexrav), and two line com- 
ponents also give fit parameters fully consistent with those mea- 
sured in the "averaged ratio flux spectrum" fitting, as shown in 
Table |2] Also in this case, a strong dependence of the broad 
line significance on the adopted continuum parametrization is 
present. 

4.3. Complex fits 

We next examine the implications of using a self-consistent re- 
flection spectrum on the broad line detection significance and its 
m easured parameters. We employed the neutral reflection model 
of iNandra et al.l (|2007|) . This model, known as pexmon, enables 
modeling of the iron lines and Compton reflection continuum in 
a self-consistent manner In addition, the model includes a neu- 
tral Fe K/3 line (7.05 keV) with a flux of 1 1.3 % of the Fe Ka fine, 
a Ni Ka with 5 % of the flux and the Fe Ka Compton shoulder 

Fitting the "averaged ratio flux spectrum" with a power law 
and a distant neutral reflector (pexmon) yielded an acceptable fit 
(X^/dof - 33.6/20). However, the fit leaves residuals around the 
Fe core indicative of the presence of a broad component. When 
this residual feature is modeled with a diskline with param- 
eters shown in Table [3] the fit improves significantly (Ax^ ~ 
1 1.0). The disk inclination is constrained to be 31.9^l^^ degrees, 
while the diskline EW is eV. We determined the reflec- 
tion fraction Roisi = 0.28 +0.06 for the reflector The top panel of 
Fig|2]shows the spectral fit (mo pow+pexmon+diskline) to the 
"averaged ratio flux spectrum". The residuals (with the pexmon 
and diskline normalizations set to zero for the sake of illustra- 
tion) are shown in the bottom panel. 

The broad residual feature was also parametrized using an 
emission line profile from an accretion disk around a maximally 
rotating black hole (laor, lLaoijll99ll) . The peak energy, Emsk 
was fixed at 6.4 keV, while the inner disk radius, Ri„ and outer 
disk radius, Ro„i, were kept at the default values of 1.24 and 
400 Rg, respectively, where Rg is the gravitational radius. The 
emissivity index, /3, was fixed to 2.5. The disk inclination to 
the observer, ;, and the line normalization were allowed to vary. 
This spectral characterization (mo pow+pexmon+laor) yielded 
an acceptable fitix^jdof - 22. 1/18) and we recovered the laor 
line EW = 112;^^!^' eV, the disk inclination / = 43.9;^g j degrees 
and the reflection fraction Roisi - 0.26 + 0.05 for the reflector 
The fit results are presented in Table |3] We note that the parame- 
ters of the illuminating power law continuum (e.g. photon index 
F and normalization) and pexmon (e.g. photon index F, cut-ofF 
energy, Fe abundance, normalization and the derivation of the 
reflection fraction Roisi etc.) were treated as in the pexrav fits, 
with the only difiference of inclination of the pexmon now fixed 
to 60 degrees. 

To model the broad feature self-consistently, we assumed 
that the accretion disk is in a low state of ionization and refit- 
ted the "averaged ratio flux spectrum" by replacing the laor 
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Fig. 7. Top panel shows spectral fit to the "averaged ratio flux 
spectrum" using a model including a power law, a distant neu- 
tral reflector and a diskline (mo pow-i-pexmon-i-diskline) along 
with the model components. Residuals (with the pexmon and 
diskline normalizations set to zero for illustration purposes) 
are shown in the bottom panel. 
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Fig. 8. Top panel shows the spectral fit to the "averaged ratio 
flux spectrum" using a power law and dual reflection spectra (mo 
pow-(-pexmon-t-kdblur*pexmon) along with the model compo- 
nents. Residuals (with the normalizations of both the reflectors 
set to zero for illustration purposes) are shown in the bottom 
panel. 



line with a blurred reflection component. We do this by con- 
volving the pexmon model with kdblur, a convolution model to 
smooth a spectrum by relativi stic effects from an accretion disk 
around a rotating black hole (iFabian et al.ll2002l) . The inner ra- 
dius of the blurred reflector (accretion disk) was fixed to 7?,„ - 
6 Rg, whereas the outer radius, Ro„t, and the emissivity index,/?, 
were fixed at 100.0 Rg and 3.0, respectively. The free parame- 
ters were the normalizations of both the reflectors (pexmon and 
kdblur *pexmon) and the inclination, /, of the blurred reflector 
The average parameters for this model are presented in Table [3] 
We find that the blurred reflector covers a solid angle {IttxRbi,,,) 
of (0.80 + 0.24):7r at the X-ray source, while the distant neutral 
reflector subtends a solid angle of (InxRoisi) of (0.54 + 0.12);Tat 
the X-ray source. The average inclination derived for the blurred 
reflection component is 3^.^^^ ^ degrees. In the top panel of Fig. 
[8]we display the spectral fit to the "averaged ratio flux spectrum" 
using a power law and the dual reflection spectra. The residuals 
(with the normalizations of both the reflectors set to zero for il- 
lustration purposes) are shown in the bottom panel. 

Considering that the inner accretion disk is likely to be ion- 
ized, we then attempted to model the broad Fe K feature of the 
"averaged ratio flux spectrum" using the blurred ionized reflec- 
tion component. For doing this, we convolved th e ionized reflec- 
tion m odel (denoted by reflionx in XSPECQ) of lRoss & FabianI 
(l200l with kdblur. The incident power law photon index of 
the reflionx model was set equal to that of the primary power 
law component and the Fe abundance was fixed to 1. The ion- 
ization parameter and the normalization of the reflionx were 
allowed to vary. The other model components such as the power 
law, pexmon and kdblur etc. were set as in the previous fits. 

This model (mo pow+pexmon+kdblur*reflionx) also 
produced a reasonably good fit with (x^/dof - 23.5/17). 
However, the f of the ionized reflector could not be constrained 
and we inferred an upper limit of 560 erg cm s ' at 68 %. 

^ In this model, an optically thick accretion disk is illuminated by 
a power law, producing fluorescence emission lines and the reflected 
continuum. The parameters of the reflionx model include the incident 
power law photon index T, Fe abundance and ionization parameter ^, 
defined as ^ = AnFj-oi/ri/j, where Fjo, is the total illuminating flux and 
nn is the density of the reflector. 



We outline the fact that the reflionx model does not include 
a formal parameter for the reflection fraction of the accretion 
disk and provides only the normalization of the reflionx com- 
ponent. Our measured value of the reflionx normalization is 
(8.2^^-^) X 10""*. We find that the parameters for the power law, 
distant reflector represented by pexmon and kdblur components 
are fully consistent with those obtained in the fit including the 
blurred neutral reflection component. The fit results are reported 
in Tabled 

In the next fits, we replaced the kdblur compo- 
nent with the kerrconv rel ativistic convolution model of 
iBrenneman & Revnoldsl (l2006h and investigated the variation in 
the fit parameters for the model comprising dual neutral reflec- 
tion components. The kerrconv model consists of seven pa- 
rameters: the emissivity indices for the inner and outer disk sep- 
arated by a break radius, inner and outer radii for the disk emis- 
sion, the spin parameter of the black hole and the inclination an- 
gle of the disk with respect to our line of sight. For the kerrconv 
component, we assumed a single emissivity index of 3. The 
black hole spin was fixed to 0, while the inclination angle of 
the disk was kept as a free parameter. We further fixed the inner 
radius of the disk to be equal to the radius of the marginally sta- 
ble orbit, Rms{= 6Rg for a non-rotating black hole), and the outer 
radius to 400 /?„,j. The best-fitting parameters for this model (mo 
pow-i-pexmon+kerrconv*pexmon) are given in Table [3] We re- 
port that the fit statistic as well as the fit parameters resulted from 
this model are in excellent agreement with those found in the fit 
using kdblur convolution model. A similar parametrization is 
obtained when the black hole spin is fixed to 0.998. 

All the above discussed complex spectral models were also 
applied to the "averaged X-ray spectrum". The fit parameters 
retrieved from the "averaged X-ray spectrum" spectral fitting 
are summarized in Table [3] It can be seen that despite the 
slightly different underlying continuum of the "averaged X-ray 
spectrum", the parameters are fully consistent with those mea- 
sured from the "averaged ratio flux spectrum". For illustration 
purposes, two spectral fits to the "averaged X-ray spectrum" 
are displayed in Fig|9] The left panel plots spectral fit using 
a model including a power law, a distant neutral reflector and 
a diskline (mo pow+pexmon+diskline), while the right panel 
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shows the fit using a power law plus dual reflection spectra (mo 
pow+pexraon+kdblur*pexmon). 

Finally, considering that the broad line EW or equivalently 
the reflection fraction of the blurred reflector, RBiun is the criti- 
cal parameter of our work and many fit parameters (e.g. cut-off 
energy of the primary power law, Fe abundance, emissivity in- 
dex, inner and outer radius of the disk etc.) were fixed in our 
spectral analysis, we investigated the effects of relaxing these 
fixed parameters on the inferred value of the RBiur- We refitted 
the data by allowing the fixed parameters to vary one by one. We 
find that the /?b/,„ parameter measured from these new spectral 
fits is fully consistent with the values reported in Table [3] and it 
is < 1.5 at the 5 sigma level. 

5. Discussion and perspectives 

It is critical to have multiple independent methods to measure the 
average Fe Ka line properties. We have determined the Fe line 
parameters using two methods described in Sect. [3] We find that 
the average narrow core EW for the whole sample comprising 
248 AGNs is ~30 eV and is stable in various continuum models 
employed for the "averaged ratio flux spectrum" and "averaged 
X-ray spectrum" (Tables[T]and|2]i. We have demonstrated that the 
broad line parameters and its detection significance are highly 
sensitive to the adopted continuum modeling. 

In order to examine the average Fe line properties as a func- 
tion of the hard X-ray luminosity (2-10 keV) and also to inves- 
tigate the variation in the line parameters obtained from the two 
methods we splitted the sample in 3 luminosity intervals (LI, 
41 < logLx < 43.5 erg s"'; L2, 43.5 < logLx < 44.5 erg s"'; 
L3, 44.5 < logLx < 46.0 erg s"'). We then computed the "aver- 
aged ratio flux spectra" from the averaged ratios and "averaged 
X-spectra" in these luminosity bins and performed their spectral 
analysis. Fig.[TO]plots the narrow Fe Ka EWs measured from the 
spectral fitting of the "averaged ratio flux spectrum" (red) and 
"averaged X-ray spectrum" (green) against the X-ray luminosity 
in the 2-10 keV band for the whole sample of 248 AGNs (hor- 
izontal lines) and its sub-samples used for luminosity bins. An 
excellent agreement between the narrow line EWs determined 
from the two independent methods is evident. The narrow line 
EW decreases as the X-ray luminosity increases, a trend des- 
ignated as the "X-ray B aldwin effect" or "Iwa s awa-T aniguchi 
effect", first reported by llwasawa & Taniguchil (Il993h . We do 
not concentrate on quantifying the X-ray luminosity dependence 
of the narrow line EW as it has been discussed in our previous 
work (see CHIO). However, we point out that the narrow line 
EWs derived in the rest-frame are consistent with the ones mea- 
sured in the observed frame. The neutral, narrow Fe Ka fluores- 
cent emission line, peaking at 6.4 keV is reported to be a com- 
mon and do minant feature in the X-ray spectra of local and dis- 
tant AGNs dYaqoo b & Padmanabhan"2004^ 'Nandr a et al.l[2007t 
IShu et al.ll2010H Corral et al. 2008; Chaudhary et al.l l201ol) . Our 
analysis further confirms this important finding. 

When the broad line parameters such as its shape, EW and 
its detection significance in the "averaged ratio flux spectrum" 
and "averaged X-ray spectrum" of the total sample of 248 AGNs 
are considered, we notice that these parameters are dependent 
on the assumed continuum and adopted stacking method as il- 
lustrated in Fig. [TT] In particular, despite having a well defined 
sample with reasonable statistics (net counts - 198 000) in the 
"averaged X-ray spectrum", we do not detect a clear extended 
red-wing, and the measured EW of the broad feature is always 
lower than 170 eV, in a good agreement with the recent stud- 
ies of the average Fe K emission from AGNs (iGuainazzi et al.l 




Fig. 10. Equivalent width of the narrow Fe Ka line measured 
from the spectral fitting (mo pow-i-gaussian-i-diskline) of the "av- 
eraged ratio flux spectrum" (red) and "averaged X-ray spectrum" 
(green) of the total sample comprising 248 AGNs (horizontal 
lines) and the sub-samples used in 3 hard X-ray luminosity bins 
illustrated by the dotted vertical lines. The EW is plotted against 
the median luminosity in each bin. Errors and upper limits are 
drawn at the 68% confidence level. 
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Fig. 11. Comparison of the equivalent width of the broad Fe Ka 
line measured from the spectral fitting of the "averaged ratio 
flux spectrum (red) and "averaged X-ray spectrum" (green) of 
the total sample comprising 248 AGNs. Respective determinis- 
tic model is listed in blue. Errors and upper limits are drawn at 
the 68% confidence level. 



120061: ICorral et al.l 120081: iLondnotti et al.ll2008h . Although we 
do not present the details of the fits in the 2-10 keV band, we 
have performed those fits. We note that the broad line EW mea- 
sured from the fits performed in the whole 2-10 keV range is 
always lower than 160 eV, which is statistically consistent with 
the value found in the 3-10 keV band. 

Our derived broad line EW is significantly lower than the 
~500 eV report ed in the stacked spectr um of AGNs in the 
Lockman Hole (IStreblvanskaetalJl2005L hereafter S05). We 
briefly discuss possible causes for the difference. Considering 
that the fluorescent line profile originating from an accretion 
disk is a strong function of the disk ionization state dFabian et al.l 
i2000) and hence of the X-ray luminosity, we first investigated 
if the observed discrepancy in the Fe K line EW is caused by 
the different luminosity distribution of the two samples. The 
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Fig. 9. Left panel shows spectral fits to the "averaged X-ray spectrum" using a model including a power law, a distant neutral 
reflector and a diskline (mo pow-i-pexmon+diskline). Right panel displays the spectral fit including a power law and dual reflector 
spectra (mo pow+pexmon-t-kdblur*pexmon). In both figures the bottom panels show respective data/model ratios. 



Table 1. Results of the "averaged ratio flux spectrum" spectral fitting in the 3-10 keV 



Model 




Power law/pexrav 


Gaussian 








Diskline 












r 


foldE Rois, 


cosi 








Eoisk 






/ 


EW 


xVdof 






(keV) (n/27T) 


(keV) 


(keV) 


(eV) 




(keV) 






(deg) 


(eV) 




Pow -1- Narrow 


1.76 H 


t 0.02 




6.40* 


0.01* 


45 ± 


13 












31.3/19 


Pow + Broad 


1.77 H 


t 0.02 




6.40* 




75 ± 


40 












25.6/18 


Pow + Disk 


1.78" 














6.40* 


6.0* 100.0* 


-2.5* 


35.3" 


128 ±41 41.6/18 


Pow + Narrow + Broad 


1.78 H 


t 0.02 




6.40* 


0.01* 


28 ± 


20 






















6.40* 




73 ± 


43 












22.8/17 


Pow + Narrow + Disk 


1.78 H 


t 0.02 




6.40* 


0.01* 


39 ± 


14 


6.40* 


6.0* 100.0* 


-2.5* 


44.6!^;° 


102!^^ 


22.1 111 


Pow -1- pexrav + Narrow -I- Broad 


1.8* 


200* 0.27 ± 


0.17 0.9* 


6.40* 


0.01* 


28 ± 


20 






















6.40* 


"•-^ -0.20 


< 105 












21.6/17 


Pow -1- pexrav -I- Narrow -I- Disk 


1.8* 


200* 0.24 ± 


0.17 0.9* 


6.40* 


0.01* 


37 ± 


14 


6.40* 


6.0* 100.0* 


-2.5* 


45.5 ± 3. 


5< 137 


21.5/17 



Notes: All en'ors and upper limits refer to the 68% confidence range for a single parameter. 

In all the spectral fits, the errors in the Fe lines EW are calculated using XSPEC eqw with err option. 

* denotes fixed parameter. 

" Error calculation is not possible due to poor fit. 



Table 2. Results of the "averaged X-ray spectrum" spectral fitting in the 3-10 keV 



Model Power law/pexrav Gaussian Diskline 

r foldE Rois, cosi Ekc (TKa EWk^ Eom R-m Ron, P i EW X^/dof 

(keV) (n/2;r) (keV) (keV) (eV) (keV) {R^) (Rg) (deg) (eV) 

Pow + Nai-row + Broad 1.64 ± 0.02 ... 6.40* 0.01* 28 ± 10 

... 6.40* 0.49!°;^2 92 ± 35 ... 20.6/17 

Pow + Nai-row + Disk 1.64 ±0.02 ... 6.40* 0.01* 37 ± 8 6.40* 6.0* 100.0* -2.5* 44.4 ± 2.4 128 ± 36 17.4/17 

Pow + pexrav + Narrow + Broad 1.67* 200* 0.29 ± 0.18 0.9* 6.40* 0.01* 29 ± 10 

... 6.40* 0.52+j;-^« 78 ± 40 ... 20.2/17 

Pow + pexrav + Narrow + Disk 1.67* 200* 0.25 ± 0.17 0.9* 6.40* 0.01* 36 ± 8 6.40* 6.0* 100.0* -2.5* 45. * 1 14 ± 41 17.4/17 



All errors and upper limits are quoted at the 68% confidence range for a single parameter. 

In all the spectral fits, the errors in the Fe lines EW are calculated using XSPEC eqw with err option. 

* denotes fixed parameter. 
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Table 3. Results of the "averaged ratio flux spectrum" and "averaged X-ray spectrum" spectral fitting in the 3-10 keV using complex 
models 



Model Power law Distant Reflector Diskline/Blurred Reflector 

r foldERoist i Eoisk Rou, P ' EW Rbiu, X'/dof 

(keV) (n/2n) (deg) (keV) (R,) (R,) (deg) (eV) (Sl/2n) 



"averaged ratio flux spectrum" 



Pow + pexmon + Disk 


1.8* 


200* 


0.28 


±0.06 


60.0* 


6.40* 


6.0* 


100.0* 


-2.5* 


Q7 Q+7.6 


OQ+50 




22.7/18 


Pow + pexmon + Laor 


1.8* 


200* 


0.26 


±0.05 


60.0* 


6.40* 


1.24* 


' 400.0* 


2.5* 


43.9!^;I 


J J 2+50 




22.1/18 


Pow + pexmon -I- kdbIur*pexmon 


1.8* 


200* 


0.27 


±0.06 


60.0* 




6.0* 


100.0* 


3.0* 


38.8!«;5 




0.40 ± 


0.12 22.9/18 


Pow -1- pexmon + kdblur*reflionx 


1.8* 


200* 


0.27 


±0.06 


60.0* 




6.0* 


400.0* 


3.0* 


39.7+5'^ 






23.5/17 


Pow + pexmon + kerrconv*pexmon 


1.8* 


200* 


0.27 


±0.06 


60.0* 




6.0* 


2400.0* 


3.0* 


43 9+2 




0.40 ± 


0.1223.1/18 



"averaged X-ray spectrum" 



Pow pexmon + Disk 1.65 ± 0.02 200* 0.24 ± 0.06 60.0* 6.40* 6.0* 100.0* -2.5* 38.5!^^ 102^^^... 17.5/17 

Pow -I- pexmon + Laor 1.65 ± 0.02 200* 0.23 ± 0.06 60.0* 6.40* 1.24* 400.0* 2.7 ± 0.6 43.9!^^ i(>l^_f^ ... 16.6/16 

Pow -I- pexmon + kdbIur*pexmon 1.67 ± 0.02 200* 0.25 ± 0.06 60.0* ... 6.0* 100.0* 3.0* 40.8 ±4.1... 0.49 ± 0.15 18.6/17 

Pow -I- pexmon + kdblurweflionx 1.67 ± 0.06 200* 0.27 ± 0.07 60.0* ... 1.24* 400.0* 3.0* 33.8!^^ 18.8/16 

Pow -I- pexmon + kerrconv*pexmon 1.66 ±0.02 200* 0.24 ± 0.05 60.0* ... 6.0* 2400.0* 3.0* 43.8!^^ ... 0.50 ± 0.16 18.3/17 

Notes: All errors and upper limits refer to the 68% confidence range for a single parameter 

In all the spectral fits, the errors in the Fe lines EW are calculated using XSPEC eqw with err option. 

In the fits including kerrconv component, the black hole spin is fixed to 0. 

* denotes fixed parameter 



luminosities of the Lockman Hole sample were taken from 
iMainieri et al l (l2002h . The 2-10 keV X-ray luminosity distribu- 
tions of 53 Type 1 and 41 Type 2 AGNs used in 805 are very 
similar, with the mean luminosity of 10^"^ *' ergs s"'. The mean 
2-10 keV X-ray luminosity of our sample is lO'*"^'^ ergs s"'. We 
conclude that the luminosity difference does not contribute to the 
difference of the line strengths. Our measured broad line EW de- 
termined from a simple power law continuum parametrization is 
much lower than the associated value derived in S05. This im- 
plies that the continuum modeling is also not the key parameter 
behind the discrepant results. The observed discrepancy is most 
Ukely due to the difference in the adopted averaging techniques. 

We also performed complex fits presented in Sect. 14.31 by 
modeling the "averaged ratio flux spectrum" and "averaged X- 
ray spectrum" with dual reflection models (a distant reflector and 
a disk reflector) handling emission lines and the Compton reflec- 
tion self-consistently. Two types of reflection models (i.e. neu- 
tral and ionized) were applied to the data. We stress that from 
the statistical view point, in both the explored averaging proce- 
dures, all models including two line components or dual reflec- 
tion spectra provide good description of the data. Interestingly, 
characterization of the distant reflector with the pexmon model 
results in a high detection significance of the observed broad fea- 
ture as compared to the pexrav. Given the available data quality 
and the well known degeneracy of accretion disk parameters, we 
could not derive all fit parameters independently. However, we 
emphasize that we have explored the parameter space in detail 
by relaxing the fixed parameters one by one. We note that the 
measured parameters from these fits are fully consistent with the 
values listed in Table |3] In particular, based on this parameter 
space exploration, we can safely say that the Rsiur parameter de- 
termined from our data is < 1 .5 at the 5 sigma level. Furthermore, 
the derived average parameters of the dual reflectors such as the 
solid angle covered by the reflectors at the X-ray source and the 
inclination etc. agree well with the average values reported in the 



local AGNs samples with high coun ting statistics dNandra et al.l 
120071: Ide La Calle Perez et a"l1l2010l) . 

Ac cording to the detailed calculations of iGeo rge & FabianI 
([T99Th . the predicted EW of the Fe Ka arising from a centrally 
illuminated cold disk is 130-140 eV for a continuum with F in 
the range 1.8-1.9 and an inclination of 40 degrees. Our observed 
mean diskline EW of ~ 100 eV agrees well with the theoretical 
prediction. Our results are also consistent with predictions of the 
relativistic Fe Ka lin e intensities from the integrated spectra of 
AGNs presented by (lBalIantvndl20 1 Ol see Figure 1). In partic- 
ular, the measured broad line EW and accretion dis k reflection 
fractio n (/?b/,„ ~ 0.5) in our data are consistent with lBallantvnd 
( 1201 Oh predictions, despite the slightly different assumptions on 
the relativistic profile for the line (laor2 versus diskline). 

In future work we plan to substantially improve the statis- 
tics by assembling larger AGNs samples. In particular, includ- 
ing MOS 1 and MOS2 data and collecting multiwavelength data 
to study the broad Fe line properties as a function of the source 
physical properties, such as the har d X-ray luminosi ty, black 
hole mass and accretion rate (see e.g. lBallantvnel20 1 Ol) . Further 
progress in exploiting the full diagnostic potential of the Fe Ka 
line in measuring black hole properties (e.g. mass and spin) re- 
quires more data with newer, more sensitive instruments with 
very high effective area and broad band coverage. The next gen- 
eration X-ray observatory, ATHENA will provide these capabil- 
ities. Owing to its huge effective area, ATHENA will enable the 
measurement and characterization of broad lines in individual 
AGNs covering a broad redshift range (z = 1-2). 
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Fig. A.l. Top panels in each figure show spectral fit to the "averaged ratio flux spectrum" using the quoted model. In the bottom 
panels ratio with respect to the quoted model is shown. The model components are also displayed. 
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Fig.A.2. Top panels in each figure show spectral fit to the "averaged X-ray spectrum" using the quoted model. In the bottom 
panels ratio with respect to the quoted model is shown. The model components are also displayed. 



